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5UMMABT 

An  analysis  of  a  solid  propellant  ramjet  which  carries  the 
propollant  in  the  combustion  chamber  was  made,  and  a  procedure 
outlined  for  calculating  the  required  burning  rate  of  such  a  ramjet 
which  develops  constant  thrust.  The  two  factors  which  influenced 
this  developnent  were: 

1)  The  solid  propellant  ramjet* s  combustion  chamber  varies 
in  sise  with  time  due  to  burning  of  the  fuel. 

2)  The  3olid  propellant  ramjet  develops  constant  thrust. 

These  factors  wore  coordinated  with  the  analysis  of  the  internal 
flow  system  of  the  ramjet,  rad  an  expression  for  the  required 
burning  rate  was  derivod. 
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I.   DPPHOUUCTIOS 

The  use  of  solid  fuel  in  b  ramjet  introduces  a  new  problem. 
As  the  fuel  bums,  the  effective  combustion  chamber  area  increases. 
This  would  normally  result  in  a  change  in  thrust?  however,  if  the 
condition  cf  constant  thrust  is  imposed,  then  this  change  in 

;tion  chamber  area  would  vary  the  amount  of  heat  required. 
The  ohaage  in  heat  required  nay  be  expressed  in  terms  of  the  burning 
rate  required  for  a  solid  propellant  ramjet  developing  constant 
thrust. 

Heferences  1  and  2,  which  are  representative  of  the  numerous 
analyses  in  this  field,  were  concerned  with  the  internal  flow  through 
a  ramjet  having  a  constant  combustion  chamber  area,  and  no  restriction 
on  the  thrust.  These  analyses  were  developed  for  subsonic  conditions, 
and  a  limited  heat  addition. 

In  the  present  Investigation,  a  theoretical  an- lysis  vr-.s   mad© 
for  supersonic  conditions  of  a  solid  fuel  ramjet,  which  carries  the 
propellant  in  the  combustion  chamber,  developing  constant  thrust. 
This  latter  requirement  fixes  the  exit  conditions  and  permits 
calculation  of  the  internal  flow  from  the  rear  of  the  ramjet  to  the 
combustion  chamber  inlet.    In  this  manner,  the  pressure  loss 
across  the  combustion  chamber  and  nos2le  was  computed.  The  pressure 
rise  across  the  diffuser  may  be  determined  from  the  Imown  inlet 
conditions.  3y  equating  the  pressure  rise  to  the  pressure  loss, 
unique  determination  of  the  heat  required  for  constant  thrust  can  be 
made.  2he  amount  of  heat  developed  is  a  function  of  the  fuel  and 


ita  distribution,  sad  tn  3,  the  required  burning  rate  nay 

)t  .'ivined. 

•  tt  of  tl  '        i  \jsl9  ar«  pretexted  i^.  Part  II  of  the 
inve  tJ     felon  in  the  form  of  vprkiag  equations,   soho  of  vhioh  Lave 

solTri  presented  on  ch-rta.     Biese  equ-'.tiona  and  charts 

hav3  oonsolldAted,  zn&  the  required  burning  r?.to  or  si>ecific 

fuel  consumption  is  i.resonted  in  Figure  10.     Part  III  illustrates 
the  method  of  solving  a  numerical  prohleE. 


-2a- 


H« 


V 


SfLCTiOt'J  A- A' 


L 


U—   D\FFU.bEft.        ->(.t-  Conflu^'ioN  CnAWitfiL 


->  A 


Figure   1    Schematic  of  a  solid  propellant  ramjet  which 
carries  its  fuel  in  the  combustion  chamber. 
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II.  AHALY3I3  I  "jOW  SYSTEMS 

Jlow  characteristics  of  the  internal  flow  system  oust  be 
calculated  so  that  the  state  of  air  entering  the  conbustion  chisabor, 
the  flow  change  across  the  combustion  chamber,  the  net  drag  of  the 
internal  flow,  and  the  required  geometry  of  the  ramjet  can  be 
determined.  In  determining  the  internal  flow  from  the  rear  end  of 
the  ramjet,  the  following  relationship  a  mast  be  established! 
(A)  nozzle  throat  condition  to  ambient  conditions  (B)  flow  changet 
across  the  nozzle;  and  (C)  flow  changes  across  the  combustion 
chamber.  Shese  expressions  must  then  be  related  to  (d)  flow  change* 
acro3s  the  diffuser. 

One  dimensional  flow,  uniform  velocity  distribution  across  the 
duct,  constant  gaa  properties,  and  Ideal  expansion  were  assumed. 

Initial  quantities  assumed  to  be  known  are  as  follows* 

(1)  She  desired  thrust  at  designed  altitude. 

(2)  The  heeting  value  and  density  of  the  fuel. 

(3)  designed  K?».ch  number  (supersonic)  and  atmospheric 
conditions  p0,  f>0,   T0. 

(4)  Skin  friction  coefficient. 

(5)  Area  of  the  throat  of  the  nozzle. 

(6)  Patio  of  length  of  chamber  to  throat  diameter. 

(7)  The  number  of  grains  (K  )  of  fuel. 


< 


A  co:.Tplote  Investigate  o  \  shoolA  Include  consideration  of: 
1)     i3entropic  diffusion,   and  2)     normal  9hock  wave;  and  whether 
the  nozzle  is  sonic  or  not  at  the  throat,     This  analysis  outlines 
in  detail  the  case  of  isentropic  diffusion  with  a  sonic  nozzle. 

.to  oicaot  relation  between  no33le  exit  condition  and  that  of 
the  anbient  surroundings  nay  he  obtained  frca  the  ©jcrression  for 
thrust.     The  thrust  of  a  ranjet  is  siren  "by  the  increase  In 
oonentga  of  the  internal  flow,  or 

y=  n  (  v4  -  vQ). 
Hie  bass  rate  of  floy  at  the  nozzle  throat  is 

ass     P\     vt     SBft     at    *S     since"tasl 
Ihe  thrust  nay  then  be  expressed  as 


or 


*  pm  Ax  lo    \  ' 


*r*i, 


_  .^[3-  m.-MoH 


-      T>. 


I'o     I     h 


•f- 


All  of  the  factors  on  the  left  hand  side  of  the  equation  are 
known  initially,  therefore  a  new  parameter,   E    ,  la  introduced 


n  by  the  perfect  gas  law 


n*    ~    p«    Art  I 


-  v-i 


ind 


»  +  I 


—     \~a-U-0 


Using  the  isentropic  relation 


y-n 


then 


E  - 


*-^< 


jr+  i 


_y-f  i 
2.  Or-') 


L  i-e^/fl;. 


AW  -  Mo  1  "^ 


L 


(1) 


Solving  for 


7; 


a- 

Ha-  -  E  L^^/V. 


M, 


(1A) 


Sj  is  not  a  known  parameter,  but  it  can  be  related  to  the  temperature 
rise.  Tli©  heat  liberated  is  equal  to  the  total  temperature  rise,  or 


l3t  -  lit 


** 


CpTzt 


(2) 


The   compression  from  free  stream  to  station  2  at  the  combustion 
chamber  inlet  is  an  isentropic  one,  and  the  state  2  is  assumed  to 
"be  a  stagnation  state,  riiere 


nt 


Similarly  the  expansion  from  state  3  to  state  4  is  an  isentropic  one 
and 


/3t   =   T^tr 


Equation  (2)  becomes 
7^r 


(at 


C^  lot 


(3) 


tiling  the  isentropic  relation 


ifct  =   7}  0  +  r^O 


equation  (3)  "becomes 


(4) 


Equating  (1A)  and  (4),  and  solving  for 
following  expression, 


_4H 


otr 


,  results  in  the 


AH 


M, 


M, 


7.»t  I 

|H«r-0 


(  +  *£** 


(5) 
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TIil3  equ-tion  c 
"be  aolvod  In  tho  form 
of  a  graph  and  would  "bo 
Binilar  to  that  shown 

In  nrnpo  2, 


\     4H 

V5.    M* 

\  CP  '•*" 

A»_ 

Cp  T»r 

Figure  2 


■The  ideal  valus  of  heat  required     (  c  -o^  )^^07  a  rriajet  yith  no 

losses  cm  be  obtained  by  entering  Hgure  2  vrtth  M4  ■  IL«     Aa  loaaes 

dovelop,  the  o:sit  mach  ntmber  H4  decreases,  md  a  larger  value  of 
AH 


Cp  ior 


is  required  to  maintain  const.snt  thrust. 


ftt,    ,  Pov  QhaPfies  Acrpua  jfoa  ff&iHfta 

In  development  of  the  analysis  from  the  rear  of  the  ramjet,  the 
flow  aero aa  the  noszl©  must  be  expressed  in  terms  of  the  combustion 
chamber  outlet.     The  pressure  ratio  across  tho  nozzle  ia 

P3 


=  f'z-  f 


T./ 


(6) 


or,  in  terms  of  mach  numbers, 


Pa 


1°. 


tin, 


r 


hence 


1+   *t*  Mz 


1+'^^ 


+  ^ 


r^- 


(7) 


1&  tho  .-.ontinulty  eqiv.tion 


fT>3    —      tnr,^ 


or 


r    /<L.  \* 


^  -  ^7— (Va_) 


-dao. 


M: 


-(*n*> 


1i?£< 


Solving  for 

<♦. 

• 

iC*-0 

T3 

■( 

^3 

If*  a-) 

r       »"»-l 

(8) 


Substitution  of  this  relation  Into  equation  (6)  yields 


i-VA* 


r-t- 1 


M|"  (/F  M 


r+/ 


(9) 


A  relation  "between     3r\kt      and  Hg  raay  be  obtained  'by  equating  (7) 
and  (8) 


^Ti>< 


N- 


*±i  \lCT-.) 


3      \     I      +     «=*    /^ 


5T  '-I-J 


(10) 


-9- 

lu  on  (9),   the  torn      y7TZ><_       is  a  known  quantity  and 

h  3  ptgrsloal  significance,  whereas  M3  13  an  unknown  value  and  is 
not  >lly  used  as  a  parameter.     This  analysis  will  therefore 

aloped  in  teito  of    ft?  x><_     I  however,  because  of  ease  of 
conversion,  as  seen  in  equation  (10),  solutions  to  all  equations 
were  derived  in  terns  of  Mg.     Die  final  result  is  expressed  in  the 
known  parameter,    fiy£<.     • 

equation  (9)  may  be  written  as 


(ii) 


Sit B.QW  Changes  Acroflg  gtt  qwtoM*toft  SW^ff 

The  flow  through  the  conbustion  chamber  is  considered  nest.     C3i© 
channel  or  duct  has  a  constant  cross  section  in  which  friction  losses 
occur.     As  the  fuel  burns  uniformly,  the  cross  section  of  the  duct 
increases,  and  a  change  in  heat  added  is  required  to  maintain  constant 
throat.     In  koepiiiQ  with  the  procedure  outlined,  all  relations  will 
be  stated  in  terms  of  combustion  chamber  outlet  conditions. 

The  pressure  drop  between  stations  3  and  3  (at  combustion 
chamber  inlet  and  outlet)  is  Obtained  by  applying  the  momentum  relation 
between  these  stations: 


P*  +  fl^   -  Pa   +-   ^  +■ 


£t-  (12) 


-2 

where  Df  Is  the  drag  force  due  to  friction  within  the  combustion 
chamber,  and  Ag  ■  Aq  . 

From  equation  (12)  it  is  evident  that  the  pressure  drop  ia 
osed  of  a  component  associated  with  tho  frictions!  drag  force, 
and  a  coriponent  due  to  the  momentum  increase  in  the  duct,  "s 
pointed  oat  in  Reference  (l),  the  former  component  is  frequently 
neglected  but  becomes  very  important  when  the  velocity  in  the 
duet  is  large.  Ihe  drag  term  is  presented  in  the  form  of  a  drag 
coefficient,  Cj>f  ,  where 


C 


**  ~  4*.  A*  (13) 


If  the  incompressible  flow  is  taken  as  a  basis  for  the  unheated  low 
speed  condition,  the  dra^  coefficient  becomes 

and  is  determined  by  the  geometry  of  the  ramjet  and  the  Reynolds  number. 
There  is  on  increase  in  the  dynamic  pressure  'rfhen  heat  is  applied, 
resulting  in  an  increase  in  the  incompressible  drag  coefficient. 
If  the  dynamic  pressure  varies  linearly  yith  the  length  of  the 
combustion  chamber,  an  average  value  may  be  assumed,  and  equation  (14) 
becomes 


C 


Dfx 


If  the  incompressible  drag  is  now  expressed  in  terms  of  the  skin 


-11- 

f rlction  coefficiont  (a  loiovn  value),    Cb^     becomes 


c 


J>fA      ~ 


or 


j_  1  (j^hlliiW'T  iVn-  i 


C^fA    =     4  X  C* 


hence 


(15) 


Since  the  heating  effect  has  bean  accounted  for  in  the  basic  drag 
relation,  the  two  drag  equations  (13)  and  (15)  may  be  equated 


or 


C«  -  ^icf  (/+%)  (16) 

With  this  definition  of  the  drag  coefficient,  equation  (12)  becones 

or 


p3    =  i +  * 


<  [tc*  ('*%)+ 0-/^0] 


-13- 


or,   In  terns  of  tho  diaiioter  ratio,  since  Dc  v  rieg, 


yri&ctit+tyyk-vti 


(17) 


In  order  to  evaluate  the  pressure  ratio,   tho  density  ratio 
must  first  "be  obtained.     This  ratio  may  he  determined  by  sinult  meous 
solution  of  the  energy  equation, 

~t  +    CPTX  I-  AH  -   -£■  +  CP73 
the  continuity  equation, 

the  perfect  gas  relation, 

and  the  aomentua  relation,   equation  (12).       The  final  result  is 


fl: 


i-      r 


i+^^^+O 


ft 


1  4-  ^  /»C 

/     T      Cp/or 


or  in  terras  of 


\pTDc 


M, 


ftti  ^    dL 


l+^/M3 


AW-['«^^*^*lt^E 


=  O 
(18) 


3olutlon  of  equation  (17)  provides    .A-      ,   for, 

I© 
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P0  R,       "     fc  (I9) 


Thus  the  pressure  dron  across  the  combustion  chamber  and  no  sale  has 
been  calculated,  and  must  equal  the  pressure  rise  across  the 
diffuser. 

Solution  of  equation  (17)  also  permits  determination  of     Mz.     , 
for  by  the  continuity  equation, 

Also, 

21     Ik  (JLiS^.\_(i  +  -^ILA  (ll^LOiL 

Corabining  the  last  two  equations,  and  simplifying, 


Mx    =  J  +  *£.'  A* 


Y-l 


^li^.r)^"      X 


(20) 


The  pressure  rise  across  the  diffuser  as  computed  from  the  rear 
of  the  ramjet  vas  plotted  vs  the  combustion  chamber  inlet  n.-ch 
number  as  calculated  in  equation  (20), (see  Figure  3). 
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Mx 


M,   vs    £ 


for  constant  E,  I"ri,  -57- 
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Figure  3 


A  plot  of  conbustion  chamber  Inlet  maeh  number  and  heat 
required  for  constant  thrust  aa  developed  in  equations  (18) 

and  (30)  is  shown  in  Figure  4. 
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(n1^ 


ah 

CpT^t 


Figure  4 
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D.       Flow  Changes  Across  The  iTlffuser 

In  an  anlysis  of  a  ramjet  having  a  designed  supersonic  raach 
number,   there  are  two  typee  of  diffusion  to  be  considered  -  namely, 
isentrordc  rad  normal  shock  wave  diffusion.     In  this  investigation, 
diffusion  was  considered  to  be  isentropic,   and  the  dlffuaer  lOO^ 
efficient.     The  pressure  rise  across  the  diffu9er  is 


T±T  = 


r-i 


H-x^    / 


(21) 


,  the  pressure  ri=se  across  the  diffuser  as  computed  from  the 
front  end  of  the  ramjet  vas  plotted  in  Figure  5. 


M, 


Mt  vs.    iL     for  constant  S,  Mq,  ^5- 


Figure  5 


P* 


p0 


Unique  determination  of  Wx.        and  ~y-         for  constant 
thrust  vas  obtained  by  superimposing  Figure  5  on  Figure  3  and 
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notixy;:  the  points  of  Intersection  a3  shown  In  Figure  6. 
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Figure  6 
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Po 


& Sat^yqAnatlpq  of,  gurnlra  flftfa 

The  condition  of  constant  thrust  has  been  satisfied  by  the 
sinultaneou3  solution  of  equations  (19),    (20),  and  (31).     With  the 

value  of      Ma.     obtained  at  the  intersection  of  the  two  curves  as 

(   4H    \ 
a  parameter,   the  v-\lue  of  he?.t  required     Vcp^t  J  was  'istemined 

After  obtaining  values 


fron  Figure  4  for  a  given  value  of     {frb^ 

•  i^      ,  a  chart 
(see  Figure  7). 


A  H  ti. 

of    -c.pN,~      for  various     y^  j>c      ,  a  chart  was  nest  constructed  of 
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Figure  7 
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l  chart  was  used  later  to  obtain  the  heat  required  for  an 
elemental  volume.     In  eonslaering  an  elemental  volune,   the 
.'it  cf  fuel  conauned  and  the  tine  of  consumption  nu.3t  "be 
culatad. 
To  aoccn^lish  this,   the  ma33  rate  cf  flow  i3  written  aa 


where 


r-t  I 


Then, 


M»-0 


The  fuel-air  ratio  nay  be  determined  fron 
f  =   h 


°* 


F  =   m 

The      A    H    was  obtained  fron  Figure  6,  and  h  is  a  given  value, 

therefore 

mF    =    rr\  f 

The  final  relation  between  the  weight  of  fuel  burned,  the  tine 
of  burning,  and  the  changes  in  diameter  nay  be  obtained  by 
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con3i Coring  a  small  elemental  volume  of  fuel,     ^he  volume  (  V    ) 
at  tine    "t      is 

V,    =    LNiK   (*M  (23) 

where    A    Dc  is  the  change  in  fuel  diaiaeter.     At  tine      *  +   4t  , 


V, 


K(t*At)=    &C  (*)   t  A  be    -    -be  (*)  +      i./VTTi>c(0 


(24) 


Or,   in  terns  of     ~^—        , 


—77—  =  dL        +    Lw-rr^^cl 


V, 


L/VMbiCt^ 


Hence,  solving  for    V,    t 


Q"£>c  (t*Ai.) 


fit    b<LJ^ 


(25) 


The  '.-/eight  of  the  fuel  consumed,  Wp  ,  may  he  obtained  from 


wF  -  ft  V, 
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!   kiae  of  fuming  of  this  el  1  volune  it 

dtb   =      "^ 

A  plot  of  the  weight  of  fuel  consumed  vs.  diameter  ratio  and 
vs,  time  of  burning  was  made  in  Figures  B  and  9. 

WF  vi    -Z=hr     $**    9  •  v«ri  L  wF   vs..  *b     -for  given    1^ 


W;: 


Wf 


T?6c 


Figure  8 


Figure  9 


■So  shoy  the  curve  in  detail,  the  specific  fuel  consumption 
minus  the  ideal  specific  fuel  consumption  wmi  plotted.     Eie  specific 
fuel  consumption  is 


s  = 


3  too    f 


3<-oo  v^p 


(26) 


For  the  ideal  case  in  which  there  are  no  losses, 


vp0Ar 


Ym    = 


I  r 


AM 


vti 


(27) 


-. 


AH 
v.lue  of    cp7oc  Tvon  Fi^ire  2  entering  with 

M*   »   W0         M  a  paraueter.     Then 

3fooo   4B  m  3. 
50    =  K  F 


where  m  is  aa  expressed  in  equation  (27). 

A  plot  was  then  constructed  of  (S-  So)  vs.  tirae  of  burning 
(see  Figure  10). 
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Figure  10 
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III.     A 

In.  a  numerical  problem  to  determine  the  burning  rate  of  .'- 
solid  fuol  ram jot  (which  carries  its  fuel  in  the  combustion 
chnabdr)  developing  constat  thrust,  tho  following  data  is  given: 


i:0   =    2 

h 

m     14,230.65  6Tu//t>. 

I 

■     25 

X      =     3 

fr 

■     140.5  1B/ft3 

°f 

■     0.003 

«"     m     1.4 

^ 

■     1  ft2 

■     50 

p0      =     4.3G2 

psia. 

*o 

e      411 .7° ft 

I 

=     1758  lb 

liquation  (5)  haa  been  solved  for  these  values,  and  is  shown  in 
Figure  11.     From  this  chart,  the  ideal  value  of  heat  required  for  a 


ramjet  with  no  losses  is  1.019.     As  the  loss  s  develop,   the  exit 

AH 

-p~«C 

required. 


xaach  number  is  decreased  and  a  larger  value  of      7f~=T,_  is 


dL  M 

The  relation  between     fc  2>t       and      ^3      as  expressed  in 

equation  (10)  is  shown  in  Figure  12. 

Solutions  of  equations  (19)  and  (30)  for  various  values  of  '^J=. 

dL 

are  shown  in  Figures  13A,  13B,  and  13C.  It  is  to  be  recalled  that 
these  solutions  were  calculated  from  the  rear  of  the  ramjet. 
Equation  (31),  calculated  from  the  front  end,  is  plotted  on  these 
same  figures,  and  the  intersection  of  these  two  curves  uniquely 
determine    ~^~  and   Mi.   for  a  constant  thrust,  solid 

propellant  ramjet. 
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The  heat  required  to  maintain  this  constant  thrust  was 
plotted  vs.  the  eonbustion  chamber  inlet  rmch  number,   and  is 
shown  in  Figures  14A,  14B,   and  14C.     Using  the  uniquely  determined 
valtfe  of   Mr   as  the  entering  parameter,   the  heat  required  was 
obtained  as  a  function  of  the  diameter  ratio    \~fiy5^  )      »   and  is 
shown  in  Figure  15. 

Using  the  process  of  an  elemental  volume  as  described  in 
Section  j:,  and  solving  equation  (25),  a  relation  between  the  wei^t 
of  the  fuel  consumed  and  the  diameter  ratio  was  determined.     This 
relation^is  shown  in  Figure  16.     The  weight  of  fuel  consumed  was 
plotted  vs.  the  burning  time  in  Figure  17.     The  slope  of  this  curve 
is  the  required  burning  rate.    Unfortunately,  for  the  values  chosen, 
the  burning  rate  is  relatively  constant.     Actually,  the  required 
burning  rate  should  decrease  with  an  increase  in  burning  time,     To 
show  the  relation  between  the  weight  of  fuel  consumed  and  the 
burning  time  in  greater  detail,  the  specific  fuel  consumption 
was  calculated  from  equation  (26)  and  the  value  of  (S  -  80)  was 
plotted  vs.  time  of  burning  and  is  shown  in  Figure  18.     This  curve 
shows  a  change  of  approximately  fifty  percent  in  (S  -  8q)  over  the 
range  of  burning  time. 
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